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Computing entanglement of an arbitrary bipartite or multipartite mixed state is in general not an easy task 
as it usually involves complex optimization. Here we show that exploiting symmetries of certain multiqudit 
mixed states, we can compute a genuine multiparty entanglement measure, the generalized geometric measure, 
for these classes of mixed states. The chosen states have different ranks and consist of an arbitrary number of 
parties. 


I. INTRODUCTION 

Characterization and quantification of quantum entangle¬ 
ment [1] lies at the heart of quantum information theory, since 
its early recognition as “spooky action at a distance” [2] in the 
Einstein-Podolsky-Rosen article [3]. Moreover, it has been 
successfully identified as a key resource in several quantum 
communication protocols including superdense coding [4], 
teleportation [5], and quantum cryptography [6]. Entangle¬ 
ment has been shown to be a necessary ingredient in studying 
quantum state tomography [7], quantum metrology [8], coop¬ 
erative quantum phenomena in many body systems like quan¬ 
tum phase transitions [9], etc. Quantification of entanglement 
is also essential for characterization of successful preparations 
of quantum states, both in two party and multiparty domains, 
in the laboratories [10]. 

The notion of entanglement is rather well-understood in the 
bipartite regime, especially for pure states [11-15]. While 
several entanglement measure can be computed for bipartite 
pure states, the situation for mixed states is difficult, and there 
are only few entanglement measures which can be computed 
efficiently. The logarithmic negativity [14] can be obtained 
for arbitrary bipartite states, while the entanglement of for¬ 
mation [12, 13] can be computed for all two-qubit states. 
The situation becomes complicated even for the pure states 
when the number of parties increase. However, there have 
been significant advances in recent times to quantify multi¬ 
partite entanglement of pure quantum states in arbitrary di¬ 
mensions [1]. They are broadly classified in two catagories — 
distance-based measures [16-19] and monogamy-based ones 
[6, 11, 20, 21]. On the other hand, quantifying entanglement 
for arbitrary multiparty mixed states is still an arduous task. 
Recently, experiments by using photon polarization [22] and 
ions [23] have been reported in which multiparty states of the 
order of ten parties have been created successfully. Such phys¬ 
ical implementations demand a general tool to compute mul¬ 
tiparty entanglement measures for arbitrary mixed states. Re¬ 
cently there have been notable advancements in this direction 
[24]. Moreover, when an entanglement measure can only be 
evaluated for pure states, the entanglement-assisted study of 
cooperative phenomena becomes restricted to only a system 
which is at zero temperature. 

We address here the question of computing the generalized 
geometric measure (GGM) [19], a genuine multiparty entan¬ 
glement quantifier, for mixed states. The GGM of pure states 
has already been computed efficiently in several systems for 
arbitrary number of parties [25]. In this paper, we define the 


GGM for mixed states via the convex roof. To deal with the 
obstacle of evaluating the convex roof extension, we use sym¬ 
metry properties of certain multiparty quantum states and sim¬ 
plify the evaluation of GGM for these classes of mixed states, 
as prescribed in Refs. [26-28] (cf. [29]). Exploiting such 
symmetries, we are able to compute the GGM of different 
paradigmatic classes of mixed states having different ranks. In 
particular, we first present the exact value of GGM for certain 
classes of rank 2 and rank 3 mixed states with arbitrary num¬ 
ber of qubits. We then compute the GGM for a specific class 
of states which is a mixture of Greenberger-Horne-Zeilinger 
(GHZ) [30] and all the Dicke states [31], having a variety of 
ranks. The common property that all these classes possesses 
is that they remain invariant under the action of same sym¬ 
metric local unitary operators on each qubit. Moreover, we 
find the GGM of a class of tripartite states of rank 4 which re¬ 
mains unaltered under different local unitaries on each party. 
Einally, we show that such symmetry properties can lead to 
an exact expression of GGM for a class of multiqudit states 
having varied ranks. 

The paper is organized in the following manner. In Sec. II, 
we review the definition and the various properties of the gen¬ 
eralized geometric measure for pure states. In section III, we 
define GGM for mixed states via the convex roof construction. 
Here, we also discuss the TerhaTVolIbrecht-Werner technique 
of exploiting the symmetry of a quantum state for simplify¬ 
ing the evaluation of a convex roof extension. The same sec¬ 
tion also contains the computation of the GGM for different 
classes of mixed states. We present a summary in Sec. IV. 


II. GENERALIZED GEOMETRIC MEASURE 

A pure state is said to be genuinely multiparty entangled if 
it is not product in any bipartition. The generalized geometric 
measure (GGM) [19] (cf. [16]) of an A^-party pure quantum 
state, \' 4 > n ), is a computable entanglement measure that can 
quantify genuine multiparty entanglement. It is defined as an 
optimized distance of the given state from the set of all states 
that are not genuinely multiparty entangled. Mathematically, 
it is given by 

£{\tljN)) = 1- Al.axili’N)), ( 1 ) 

where Aj„ax(|V'Af)) = iiieix |(x|'!AAr)|, with the maximization 
being over all |x) that are not genuinely multiparty entangled. 
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An equivalent form of the above equation is [19] 

= l-max{Aj.^|/UL = {Ai,= 0}, 

( 2 ) 

where Xi-,l is the maximal Schmidt coefficient in the bipartite 
split I : L of IV’ at)- 

Let us enumerate some properties of the GGM which estab¬ 
lish it as a bona fide measure of genuine multiparty entangle¬ 
ment [19]. fdV'tv)) > 0, for all IV-jv), f(IV’tv)) = 0 iff IV’iv) 
is not genuinely multiparty entangled, and £’(|V’Ar)) is nonin¬ 
creasing under local quantum operations at the N parties and 
classical communication between them. 


III. GGM FOR MIXED STATES 

We can now define the GGM of a general mixed quantum 
state, in terms of the convex roof construction. For an arbitrary 
W-party mixed state, pM, the GGM can be defined as 

Q{pn)= min VK^dV’w))) (3) 

{pi.lV-L)} i 

where the minimization is over all pure state decompositions 
of Pm i.e., pN = Fd'*/’Ar)(V’ 7 vl' i® difficult to find the 

optimal decomposition and the computation of GGM is in 
general impossible even for moderate-sized systems. How¬ 
ever, the situation is different if the mixed quantum state un¬ 
der consideration possesses some symmetry [17, 27-29]. In 
Ref. [27], Vollbrecht and Werner have provided a general 
method to compute an entanglement measure, defined via the 
convex roof extension, of a class of mixed states which are 
invariant, on average, under a group of local unitaries. Be¬ 
low we briefly outline the same. Suppose p'm = {Ui®U 2 ® 
... 0 Um)pn{U\ ® Ul ® ® Ujf), where Ui are the local 

unitary operators, acting on Hilbert spaces Hi. The GGM of 
Pm and are the same. If it happens that pM = p'm^ then 
{Ui 0 1/2 (g)... (g) {7jy) is called a local symmetry of pM- Let G 
be a group of unitary operators U = (C/i0C/20...0f7jv) and 

P be a twirl operator, such that, A f dU UAW = P(^), 
where the integral is carried out Haar uniformly. In case of 
a mixed state pM, if there exist a twirl operator P such that 
P(pAr) = Pm, then the entanglement, G{pn), can be obtained 
from a pure \ip) which satisfies 

Pd^)(V’|) = p. (4) 

In principle, one can have a set of pure states, {IV')} = 
which satisfies Eq. (4), and it is sufficient to perform the op¬ 
timization over this set. A further step is needed where we 
convexify the optimized quantity over the parameters in pM, 
if it is not already convex. 

We now show that this method can be utilized to evaluate 
the GGM for several classes of multiparty states with arbi¬ 
trary number of parties having certain symmetries. We present 
these classes according to their ranks. 



X 


FIG. 1. (Color online.) GGM of Pm(x) = 2;|'i/’iv}(V’iv| + (1 ~ 
a;)|V’iv)(V’ivl against x. All the quantities are dimensionless. 


A. Classes of rank 2 multiqubit states 

The rank 2 mixed state, which we are now going to consider 
is a mixture of two orthogonal W-party pure states, given by 

p%{x) = x\tpM){tpN\ + (1 - x)\lpM)('ipM\, (5) 

where the subscript and superscript of p represent the number 
of qubits and rank respectively. Here, IV’Af) |'i/’)v) li® 
two orthogonal mutually complementary subspaces of the N- 

party Hilbert space . IV’Af) = with 

(^) 



i=l N-k k 


where \Dg)’s are the generalized Dicke states [31] with k 
number of excitations i.e. they are the general superpositions 
of pure states with all permutations of {N — k) |0)’s and k 
|l)’s. And 

[Nj-i 

|V’)^)= ^ a[\Df+^). (7) 

i=0 

We have chosen the coefficients in all pure and mixed states 
such that there are properly normalized. 

For pjii{x), we can find a group of local unitary op¬ 
erators consisting of two unitaries, U\ = I, and U 2 = 
(Jz, which, on average, keep p%{x) invariant. Here, I 
is the identity operator on the qubit Hilbert space and 
ax, ay, and az are the Pauli operators. One can check 
that p%{x) = ELi where 

l'0Ar(a’)) = y/xl-ipM) + — a;|V')v) is ths only class of 

pure states that is twirled to pff{x) by applying the twirl oper¬ 
ator corresponding to those unitaries. Hence, by following the 
recipe in [27], we can calculate the GGM of p%{x). Since it 
involves several parameters, for illustration, we choose fully 
symmetric states, i.e, when all the coefficients of IV’at) and 
IV’tv) are equal. The GGM of p%{x, sym) is the convex hull 
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FIG. 2. (Color online.) A plot of the GGM of p%{x\,X 2 ) = 
XI \GHZ+){GHZ+\ + X 2 \D^){D^\ + {1 - xi - X 2 )\D^){D^\ 
with the state parameters a;i and X 2 - All the axes are dimensionless. 

of the GGM of the pure states sym)) = ^/x\'4>N) + 

VT — Here the phase, gives the mini¬ 

mum GGM among all the GGM with different (j) values. We 
then hnd that GGM reaches its minimum for (pmin = 0. 
Therefore, the GGM of p\{x, sym) is given by 

Q{p%{x,sym)) = ^(1 - 2 y/x^/l - x), ( 8 ) 

since the right hand side is already convex as depicted in 
Fig. 1. An important point to note here that the GGM of 
p%{x, sym), given in Eq. ( 8 ), is independent of number of 
parties, N. 

B. Classes of rank 3 multiqubit states 

We now calculate the GGM for different classes of mixed 
states, of rank 3. 

1. Case 1 

Let us now consider a three-qubit rank 3 mixed state, 
p^{xi,X 2 ) [17], which is a mixture of known \GHZ)^), \D^), 
and |£)2). Here, \GHZ+) = ^(|000) -f |111)) [30], and 

\D^) and |D^) are given by \Dg) and \Dg) ofEq. ( 6 ) respec¬ 
tively, with bkj = for all j. It reads as 

pUxi,X 2 )=xi \GHZ+){GHZ+\+X 2 \D^){D^\ 

+ il-xi-X 2 )\D^){D% (9) 

Note that \D^) is the well-known W-state [32]. The mix¬ 
ture is invariant under local unitaries given by Ui = I, 

U 2 = 2 ri V and 1/3 = 1 „ V when they act on 

\0 e 3 y \^0 e~3— J 

each qubit [17]. The corresponding pure state which after lo¬ 
cal unitary transformations, leads to pg (xi, X 2 ), can be written 
as 

|' 0 |(xi,X 2 )) = yG^lGHZ+) + 

-FVI -xi- X 2 e^‘^^ IB^).(IO) 



FIG. 3. (Color online.) Plot corresponds to GGM of vs. the 

mixing parameters xi and X 2 - Here, a = 0.55 for the \gGHZ 3 ) 
state. Both convex and nonconvex regions are seen. The convex part 
corresponds to the GGM of p^’^{xi, 2 : 2 ). All quantities are dimen¬ 
sionless. 

The minimum of GGM among {(/)i,(/)2} is again obtained 
when 4 >i = (l>2 = 0 . By computing the Hessian matrix, 
we hnd both analytically and numerically that the GGM of 
|'03(xi,X2)) is convex with respect to Xi and X2- Therefore, 
the GGM of P3(xi, X2) is given by 

g {pI{xi,X2)) = i ^3 - |l - 5xi - 12 x2(x2 - 1) + 

8-\/6xiX2 ^1 -F sj X2(1 — Xi — X2) — xi — X2^ -F 

4 xi (^1 -F 3\/x2(1 - Xi - X2) - 8x2) | ( 11 ) 

and is depicted in Fig. 2 . 


2. Case 2 

Let us now move to a more general state while keeping the 
rank hxed. Precisely, we consider a class of mixed states of 
the form 

P3’®(xi,X2) = xi\gGHZ'i){gGHZ3\ +X2\DI){DI\ 

+ {l-X3-X2)\Dl){Dl\, ( 12 ) 

where \gGHZ 3 ) = a|000)-F'\/l — a^jlll) is the generalized 
Greenberger-Horne-Zeilinger state with 0 < a < 1. The set 
of local unitaries that keep p^xi, X 2 ) invariant, also keep the 
state P 3 ’®(xi, X 2 ) invariant, and the class of pure state that are 
projected to p 3 ’®(xi,X 2 ) is given by 

|'03’®(xi,X2)) = yGGlgGHZs) -F 

+e*^Vl-a;i -X2|Z7^). (13) 

In this case, we have p 3 ’®(xi,X 2 ) 

Ej=l Uf^\i’3^{xi,X2)){'tp3^{xi,X2)\U^g®^, 

{Uj,j = 1, 2, 3} is the same as in Case 1. 


where 












4 



FIG. 4. (Color.) Plot corresponds to GGM of IV' 3 ’®) vs. xi, for two 
values of r = Here, a = 0.55 for the \gGHZ 3 ) state. These 

are given by the dotted lines. The straight lines corresponds to the 
convexified quantities. All quantities are dimensionless. 


Numerical simulation guarantees that the minimum of 
^(l' 03 ’^(^ii 2 : 2 ))) occurs for (j)i = 4>2 = 0. However, un¬ 
like the previous cases, we find that {xi,X 2 ))) is not 

convex for all values of xi and X 2 - In particular, we plot 
^ 2 ))) in Fig- 3, when a = 0.55 and when the 
coefficients in \Dg) and \Dg) are all equal. For certain re¬ 
gions of the parameter space, the figure is already convex, 
and hence the GGM of (xi, X2))) in that region is the 
GGM of p^’^{xi,X 2 )- On the other hand, for the remain¬ 
ing regions, a convexification has to be carried out to obtain 
the GGM of p^{xi,X2)- Specifically, {xi, X2))) 

Q{p^'^{xi,X 2 )), when xi is high while X 2 is low. To obtain 
the GGM in that region, the convexification is required. To 
illustrate the process, we introduce a new variable, r = 
and let us consider cases where r = 0.96 and 0.98. The con¬ 
vexification of the curves so generated are depicted in Fig. 4. 


FIG. 5. (Color online.) The plot of GGM for p 5 (a;i,a; 2 ) = 
XI \GHZ+){GHZ+\ -f a ;2 \D^){D^\ -f (1 - si - a; 2 )|I>^)(I>^| 
against a;i and X 2 whenever it is convex. All axes are dimensionless. 

is clear that for example the GGM of a; 2 )) is convex 

for 0.64 < xi < 1.0 and 0.0 < X 2 < 0.36 and hence in that 
region, we have the GGM of pg(a;i,a; 2 ). In the rest of the 
region, to obtain the GGM of pg(a;i, a; 2 ), we have to find the 
convex hull of £{\'ij: 1 {xi,X 2 )))- 


C. Higher rank multiqubit states 

We now consider classes of mixed states with rank more 
than three. First, we explore a class of multiparty states which 
can be dealt with symmetric unitaries. In other words, this 
class of states remain invariant, when the same unitary acts on 
all the parties, i.e. p^ = Yhj We will then 

find another class of states for which symmetric unitaries do 
not work. 


1. Symmetric unitary case 


3. Case 3 


Let us move to a class of states which is a multiqubit gen¬ 
eralization of ^ 3 ( 0 : 1 , X 2 ). It is given by 

p\{xi,X2)=x^ \GHZ+){GHZ+\+X2 \D^){D^\ 

+{^-x^-X 2 )\D^-^){D^-\{U) 


where |GiTZ+) = ^(|0)®^-F |1)®^), and is given 


by \D‘-^ 1) of Eq. ( 6 ) with bkj = Again, we have 

y\k) 

P%{XI,X2) = Yl]=lUf^\tlt%{xi,X2)){tlt%{xl,X2)\U]®^ , 
where |^^(xi,X 2 )) is given in Eq. (10) with \D'^) being re¬ 
placed by for the same set of unitaries, given in Case 

1. Hence, we can compute the GGM of X 2 )) and 

check its convexity. Eor (pi = (j )2 = 0 which gives the lowest 
GGM, Eig. 5 shows the GGM of \tp^{xi,X 2 )) with respect to 
the parameters, xi and X 2 with N = 5. Erom the figure, it 


Let us now consider a class of mixed states with arbitrary 
number of parties, which can be obtained by generalizing 
pI{xi,X 2 )- The state, p^{xi,X 2 , ■ ■ ■, xn-i), is a mixture of 
generalized GHZ and all the Dicke states. It reads as 

p%ixi,X2, ■ ■ .,XN-i) = (1 -^x^)\gGHZ]sr){gGHZ]sr\ 

i 

N-1 

+ J2x,\Dl){Dl\,(15) 

with \gGHZN) = -F \/l — a^jl)®^. Rank of the 

above state spans the integers in [1, iV]. One can check that 


N 


p%{xi,.. .,xn-i) = ^ Uf^p^ixi,.. .,xn-i)U. 






( 16 ) 
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FIG. 6. (Color online.) GGM of pf = xiP\GHZ+] + ^{P[D^] + 
P[D^]) + p|£)4j^ dimensionless. 


FIG. 7. (Color online.) Plot of GGM of p| with respect to the pa¬ 
rameters, X and y. The GGM of the corresponding unique pure state, 

It/'K®, y)) = y®ICi) - ‘i^/W^{\C, 2 ) - ICs)) + Vl-*-y|C4) has 
a kink along the lines shown on the surface, in the plot. The GGM of 
the pure state is non-convex around these lines, and hence convexifi- 
cations are required thereat. 




where the set of local unitaries, consists of I and 

^ with j = 1,..., {N — 1). We have to now show 

that 

p%{xi,X 2 ,.. .,XN-l) = Uf^\il;^{xi,. . .,XN-l)) 

j 

(ij^ixu...,XN-l\Uj®^, ( 17 ) 

where \tlj^{xi,... ,xn-i)) = - J 2 i Xi\gGHZN) + 

To prove this, the we note the actions of lo¬ 
cal unitaries on each off-diagonal terms which e.g. are given 
by 


Uf^\Dl){Dl\u]^^ = e'-^^\DI){Dl\. (18) 

_ 27 ri(q — r) ^ ___ 

We use the identity e ^ = Oqr in the analysis. Simi- 

larly, 

Y^Uf^\Dl){gGHZN\U]<^^ = e^\D^^{gGUZ^\ = 0 . 

(19) 

All off-diagonal terms therefore vanish. We can now cal¬ 
culate the GGM of \ipf^{xi,..., xn-i)) and check whether 
£{\ip^{xi,..., xn-i))) is convex or not. If it is convex, then 
£(\i^§lxi,... ,XN-i))) = Gip^ixi,. ■ ■ ,XN-i))- Other¬ 
wise, we have to perform convexification to obtain the exact 
value of G{p%{xi, ... ,xn-i))- To illustrate this example, we 
consider a five-qubit state which is of the form 

pI = x^\GHZ+){GHZ+\ + ^i\D^){D^\ + \D^){D^\) 

+ l-a.l-at2 (|^3)^^3| + |^4)p4|)^ 

( 20 ) 

Following the aforementioned prescription, we compute 


£{\iI^^{xi,X 2 ))) with 


X 2 )) = ^,\GHZ+) + Jy E 1 ^') 




1 — XI — X 2 


I (21) 


k—3 


For (pk = 0, fc = 1,..., 4 which gives the infimum of GGM, 
£{\il)l{xi,X 2 ))) is plotted with xi and X 2 in Fig. 6. By using 
the Hessian technique, we find that it is convex for the entire 
range of xi and X 2 - Therefore, (/(pi) is obtained for all xi 
and X 2 and is given by 


G{pI) 


1 /1 G ^f2xi+4:X2+5 7 — 2x1 — 4x2 

2 I ~ ~ I 10 10 


/ /xiX 2 /xi(l — Xl— X 2 ) 2x2 2(1 — Xi — X 2 ) 

-^1 - 3 : 2 )) I) y (22) 


Comparing Figs. 5 and 6 with the situations obtained before, 
it seems that higher rank states, for a fixed total number of 
qubits of the entire systems, have a greater affinity for being 
convex, when their GGMs are considered. 


2. Asymmetric unitary case 

Until now, we have considered the states which remain un¬ 
altered under local symmetric unitaries of the form Uf . 

Let us now illustrate a class of three-qubit mixed states 
which remains unchanged under the local unitaries of the form 
Ui ® Uj 0 Uk- The class of mixed state having rank 4, reads 


pt = 

i 


(23) 
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where 


|Ci) = i(|001) + |010)-|100) + |lll)), 

IC 2 ) = ^HIOOO) - i|oii) +1100) + 

IC3) = i(*|000) + z|011) + |100) + |lll)), 

and IC 4 ) = i(|001) + 1010) + |100) - 
It is invariant under {Ui,i = 1,... 4}, which are given by 


Ui=I®I®I, 

U 2 = iCTy 0 H' (§1 H\ 

U^=I®ay® ay, 
and 1/4 = —lay 0 H'^ 0 H'^, 


"O'H' = ^ 1-1 1 


Note that these unitaries form a 


closed group. The only pure states that are twirled to the 


above mixed states are of the form 


XiC 




\Q- 


We compute the GGM of and minimize it over (j)i's. The 
GGM of P 4 is given by the minimum of the £’(|'0|)) for dif¬ 
ferent values of ^^s provided the quantity is convex itself. 

To visualize its GGM, let us consider, X 2 = X 3 = |, i.e. 
the state is of the form 

p1 = ^ICi)(CiI + !(IC2 ) (C2| + |C3)(C3|) 


+ (1 - X - 2 /) IC 4 ) (C4|. (24) 


In this case, we find that the minimum GGM of y)) for 

different values of is obtained when 4 >i = — (t >2 = —^ 
and 6 -i = 0. We find the GGM of pi(x,y) by convexifying 
the GGM of IV’Kx,?/))). 


D. Cases of multiqudit states 

In the previous sections, we have evaluated the GGM of 
certain multiqubit systems. We will now show that a similar 
method can be extended to obtained the analytical expression 
of GGM of multiqudit mixed states. Specifically, we consider 
an iV-qudit mixed state of rank d, in the Hilbert space , 
of the form 

d 

= (25) 

k^l 

where |T')fc = E{ 4 } 9 ^ 112 .|jij 2 -jAr)^'^^ and 

(^^jm)(mod d) = k. Our aim is to evaluate the 
GGM of the state Therefore, like previous cases, we 

construct a twirling operator, consisting of unitary operators 
Zd which are d-dimensional, non-hermitian generalization of 
the and given by 

d-l 

Zd = '^e^\j){j\. (26) 

1=0 


Here, each of the unitary operators act locally and 
symmetrically on p% d Zf^. Note that the set 

|/d, Zf^, {zf^^ ,.., {zf^^ I forms a group and the 

corresponding twirling operator keeps ^ invariant. Now, 
we have to find the set of all pure states that 

are projected to pf^ ^ under the action of the aforemen¬ 
tioned twirling operator. It can be easily checked that 
1 ^^)% d — J 2 k=i are the only class of pure states 

that are mapped to under the twirling operator, i.e., 

EtZo {zf^y\'^)%,d{'i’\%,d{zl^''y = p%,d- In this case 
also, the minimum of the GGM’s of |'I')^ ^ over the phases 
{(t>k} gives the GGM of p^ ^ provided the minimum GGM is 
already a convex function of the state parameters. Otherwise 
one has to convexify the function to obtain the GGM of p^ 
Until now, we have considered systems with the same di¬ 
mensions of the local Hilbert spaces. However, this formal¬ 
ism can be further extended where the local Hilbert spaces’ 
dimensions are not equal, i.e., for quantum systems belonging 
'^n'Hdx®did 2 ® ■ ■ ■'Zi'HdN, with di ^ d 2 ^ ...d^. In that case, 
we have two different scenarios. Firstly, oidi = 02^2 = = 

div, where . Without loss of generality, d^r is 

taken to be the largest dimension and the corresponding uni¬ 
taries are of the form ® Zd 2 -- 0 Zdfj with its subsequent 
powers upto dA? — 1, such that the composite unitary matri¬ 
ces form a group. Evidently, the case of equal dimensions is 
a special case of this. Thus, the pure states over which we 
have to perform the minimization still have the same form, 
with a slightly different version of the condition given by 
jm(inod d^r) = k. The second one is the situation when 
all the dimensions are prime to each other, and in this case, 
we have to take unitaries upto the power of (did 2 ...dAr) — 1, 
where the form of pure states remain the same, with the mod¬ 
ified condition, (mod did 2 ...dAr) = k. Therefore, in 

general, we have to take the maximum power of the unitaries 
which is the lowest common multiple of di, d 2 ,..., d^r to ap¬ 
ply the similar prescription. In the next paragraph, we illus¬ 
trate this with an example. 

For simplicity, we consider the following three-qutrit 
state, pI^ = J 2 l=oXk\^}k{^\k, where |T')fc = 

Ej and ji + j 2 + j 3 (mod 3) = k. The 

exact form of the pure states reads as 

l'I^o) = i(EN*)+ E 1012)), 

i—0 perm 

I'I'i) = E 1001) + E 1022) + E |112))> 

perm perm perm 

and 14 / 2 ) = E 1011) E 1002) + E |122))- (27) 

perm perm perm 

For this case, the unitaries which construct the twirling op¬ 
erators are given as {I^,Z^,Z^'\. Note that the unitaries of 
the form Z^ 0 Z^ Z^ form a group for i ranging from 0 
to 2 and pg 3 is evidently invariant under the corresponding 
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FIG. 8 . (Color online.) Plot of GGM of p 3_3 against xi and X 2 - The 
GGM of the corresponding unique pure state, IT's, 3 ) = y/xi\'^)i + 
— x\ — a: 2|^')3 is plotted with (j)i = (f )2 = 
0. The GGM of the pure state is convex everywhere, as evident from 
this plot and hence f ( 1 ^ 3 , 3 }) = Gip^a)- 

twirling operator. The pure state that is mapped to pg 3 un¬ 
der the action of the aforesaid twirling operator is of the form 
l^i. 3 ) = y/^l«')i+e*^=V^|4')2 + e*^Vl - a;i - a:2|«')3. 
It can be easily found that minimum GGM of [thg 3 ) is ob¬ 
tained for (j)2 = (l>3 = 0 it is a convex function of the 
parameters xi and X 2 - Hence, the GGM of pg 3 is given by 

Gipl.s) =§{1 - \/xi{l -xi- X 2 } 

- ^/X2{1 - Xi- 0 : 2 }}. 

G{p 3 3 ) is depicted in Fig. 8 and the convexity of the function 
can be visualized from the same. 


IV. CONCLUSION 


Computing entanglement of an arbitrary mixed state is a 
formidable task. The entanglement of mixed states is gen¬ 
erally dehned by constructing the convex roof over all pos¬ 
sible pure states which is practically impossible to compute 
in most of the cases. Although there exists a few bipartite 
measures which can be obtained for arbitrary states, the eval¬ 
uation of entanglement for a mixed state in multiparty domain 
is still a challenging task. In this paper, we have computed a 
genuine multiparty entanglement measure known as general¬ 
ized geometric measure of some classes of mixed states with 
arbitrary number of parties and dimensions by using certain 
symmetries. We evaluate the measure for several classes of 
multiqubit and multiqudit states having different ranks. The 
method, we exploited, uses a pure state that contains the same 
amount of entanglement as the given mixed state, and leads to 
the mixed state by action of a certain twirling operation. 


Note added: The present work is based on a poster pre¬ 
sentation [33] at the International Workshop on Quantum In¬ 
formation (IWQT2012), Harish-Chandra Research Institute, 
Allahabad, India. We thank J. Solomon Ivan for pointing out 
during a discussion over the poster that the same method as 
followed here can be used to evaluate the GGM for an arbi¬ 
trary mixture of \GHZ^) and \GHZ^), where \GHZ^) = 
-^(lO)®-^ ± |1)®^). We thank Otfried Giihne for informing 
us about their independent work on evaluating multipartite en¬ 
tanglement [34], by a method that is different from the one 
followed in the present work. 
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